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VERY ACCURATE APPROXIMATIONS FOR THE ELLIPTIC 
INTEGRALS OF THE SECOND KIND IN TERMS OF 
STOLARSKY MEANS 

ZHEN-HANG YANG 
Dedicated to the my mother Ru- Yi Jiang 


Abstract. For a, & > 0 with a ^ b, the Stolarsky means are defined by 
/ oiaP - bP) \ 

iipq{p-q)^0 

and Sp^q (a, b) is defined as its limits at p = 0 or q = 0 or p = q if pq (p — q) = 0. 
The complete elliptic integrals of the second kind E is defined on (0, 1) by 


Eir) 



— sin^ tdt. 


We prove that the functions 

1 ~ 'S'll/4,7/4 (li ^0 


and G (r) 


l-(2/7r)E(r) 

1 ~ 'S'5/2,2 (!> ''"0 


are strictly decreasing and increasing on (0,1), respectively, where r' = \/l — r^. 
These yield some very accurate approximations for the complete elliptic inte¬ 
grals of the second kind, which greatly improve some known results. 


1. Introduction 


For r G (0,1), the well-known complete elliptic integrals of the first and second 
kinds m, m are defined by 

(1.1) K{r)= I =, A (0+) = A(l-)=oo 

Jo yl — sin t ^ 

and 

/•^/2 ! - TT 

(1.2) A (r) = / V 1 — sin^ tdt, A (O'*') = —, A (l“) = 1, 

Jo 2 

respectively. These integrals can be expressed exactly in terms of Gaussian hyper¬ 
geometric function 


2 Fi(a, 6 ;c; z) = 


^ (a)„ (^)„ (a3)„ 


n=0 


(c). 


,1^1 < 1, 


Date: May 10, 2015. 

2010 Mathematics Subject Classification. Primary 33E05, 26D15; Secondary 26E60, 26A48. 
Key words and phrases. Complete elliptic integrals for the second kind, Stolarsky means, 
monotonicity, inequality. 

This paper is in final form and no version of it will be submitted for publication elsewhere. 


1 












2 


ZHEN-HANG YANG 


where a, 6, c S K with c ^ 0, —1, —2 ,(a)„ is defined by (a)Q = 1 for a ^ 0 and 
(a)^ = a (a + 1) ■ ■ ■ (a + n - 1), a 0. 


Indeed, we have 

oo /" 1 \ 2 

1 V 2n 

2 ^ n!n! 

n—0 



(1.3) 


K = 


K(r) = -F 


1 1 


2 ’ 2 ’ 




It is well known that the complete elliptic integrals have many important appli¬ 
cations in physics, engineering, geometric function theory, quasiconformal analysis, 
theory of mean values, number theory and other related fields m, a, m, i, 

0, a, m, [H], m, [la. 

Let 1(1, r) be the arc length of an ellipse with semiaxis I and r G (0,1). Then 

I (l,r) = 4E{r '), 


where and in what follows r' = '/I — r'^. In 1883, Muir [14] presented a simple 
approximation for ^ (l,r) by 2 ttA 3/2 (!,?'), where 


/nPMbP\^^^ , _ 

Ap (a, b) = I - - — j if p 0 and (a, b) = 


is the classical power mean of positive numbers a and b. This mean contains some 
simple ones, such as A-i(a,b) = H (a,b) -harmonic mean, Ai(a,b) = A(a,b) 
-arithmetic mean, Aq (a, b) = G (a, b) -geometric mean, and A 2 (a, b) = S (a, b) 
-root-square mean, etc. 

In 1996, Vuorinen |15] conjectured that the following inequality 

(1-5) E{r)>^A,/2(l,r')=--lA,(r') 

holds for 0 < r < 1. This conjecture was proved in 1997 in m Theorem 2] by Qiu 
and Shen (see also 0 Theorem 1.1]). Barnard et al. 0 discovered an upper bound 
(7r/2) A 2 (1, r') for E (r), that is, 

(1.6) E{r) <^A 2 il,r') (0 < r < 1). 

An improvement of (11.61) was presented by Qiu in (TT] Corollary (1)] (see also [TSl 
Theorem 22.]), which states that the inequality 

(1-7) E(r)<^A,„(l,r') 

is valid for all r G (0,1) with the best constant go = In2/In (7r/2). 

Motivated by the inequalities (fra-ifTTi). some new approximations for E (r) 
in terms of bivariate means were presented. For example, Chu and Wang m 
Corollary 3.2] proved the inequality 

(1-8) E(r)<^Cy,(l,r')=:^Ae{r'), 

holds for r G (0,1) with the best constant 1/4, where Cp (a, b) is the Lehmer mean 
of positive a and b defined by 

flP+i fjP+i 


Fp (u, 6) 


aP + bP 
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Soon afterwards, Chu et al. [201 Theorem 3.1] gave a much better upper bound for 
E (r). They proved the double inequality 

16-3^^ A (1 /) _ 3^-8 S ( 1 , /) 

(1.9) < -E (r) <1^(1, /) - (1, r') + (1, r') =: A, (r') 

TT 8 16 16 

holds for T e ( 0 , 1 ). 

Wang et al. [211 Theorem 2.4] established a sharp double inequality for E (r): 
OQ c; 1 o 

(1.10) A-, (r') := -Al (1, /) - -H (1, r') - g 5 (1, r') < -E (r) 

24 - 5 ^/ 2 ^ 3(1 /)_ 8 -^-%/^ ff (1 /) _ ( 38 - 3 ^) 5 ^( 1 ,/), 

27r(3-2A/2) ^ ’ ’ 27r(3-2A/2) ^ ’ 2T^(Z-2^f2) ^ ^ ’ 

and pointed out that the lower and upper bounds in fll.lOp are stronger than ones in 
and (II. 6 L respectively. In 2013, Wang and Chu HH Corollary 3.1] presented 
another improvement of dni) and (fra . which states that 

1 ("Oj-'Z I 14 / I q\2 „ ^_ 

(1.11) Asir')--—^ -— - 3 —<-E (r) < - V4r'2 + - 8) r'+ 4. 

128 (r' + 1 ) tt tt 

Very recently, Hua and Qi showed in [23l Theorem 1.3.] that the double inequal¬ 
ity 

( 1 . 12 ) 

8 


. , 1 -I- r' -I- 1 -I- r' 2 

3‘<’'>^=T(TT7r+— 


- 2 


1 + r' + r'2 / 6 


(l-fr') 


is valid for r G (0,1). 

Other various approximations for E (r) can be found in uni, ED, m, i, m, 
ESj, mi, ES], ESI; EOj: ED: E2]: E2]: and references therein. 

For a, 6 > 0 with a b, the Stolarsky means Sp^q{a, b) are defined in [5^ by 

jP -hP \ Gp 

.p(1„TTmi) itP7^o.« = o, 

(1.13) Sp^q{a,b)={ ( 

if P = 0,9 7^0, 
if p = g ^ 0, 


g(lna — In b) 

aP In a — &P In & 1 


exp 

. Vab 


aP - bP 


P 


if p = g = 0, 


and, Sp^q{a,a) = a. This family of means contains many famous means, for ex¬ 
ample, S'i,o(a,6) = L{a,b) -the logarithmic mean, S'i4(o,6) = I{a,b) -the iden- 
tric (exponential) mean, S'2,i(a, &) = A(a,b) -arithmetic mean, <S'3/2,1/2(0, &) = 
He{a, b) -Heronian mean, S'2p,p(a, b) = A^/P{aP, bP) = Ap -the p-order power mean, 
5'3p/2,p/2(0) ^) = He^f^(aP = Hs.p -the p-order Heronian mean, Spfi{a,b) = 
L^/P{A‘A^) = Lp -the p-order logarithmic mean, Sp^p{aA) = I^/P{aP,bP) = Ip 
-the p-order identric (exponential) mean, etc. 

Stolarsky means have many well properties, which can follow directly from the 
defining formula (11.1311 and be found in [35] , [3S] , [37] , [38] , [39] , [40] , [41] , [42] , [43] , 
[44] . [45]. For later use, we mentioned the following properties: 
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(PI) For all a,b > 0 and p, q S R, Sp^q{a, b) are increasing with both p and q, or 
with both a and b (see [36l (2.22), (3.12)]). 

(P2) For fixed c > 0, Sp^ 2 c-p (a, b) is increasing in p on (—oo, c] and decreasing 
on [c, oo) (see (3.14)], [431 Corollary 1.1]). 

(P3) For fixed c > 0 and p G (0,2c), {l/9p) Sp^ 2 c-p{a,b) is decreasing in p on 
( 0 , c) and increasing on (c, 2 c), where Op is defined by 


(1.14) 


Op = f- j if p ^ c and Oc = 

(see HH Corollary 1.2]). 


Remark 1. Taking a = b in (PS) yields that (1/Op) Sp^ 2 c-p {a,a) is decreasing in 
p on (0,c) and increasing on (c, 2c). Since Sp^ 2 c-p {a,a) = a, it follows that Op 
strictly increasing in p on ( 0 , c) and decreasing on (c, 2 c). 


Now we intend to estimate for the complete elliptic integrals of the second kind 
E (r) by the Stolarsky means of 1 and r', i.e. Sp^q (1, r'). Expanding in power series 
gives 





P + 9 - 9/24 p + g-9/2 g 

- r - r 

96 128 

8{p + q) (2p2 + 550 ) + 19 845 g . 4 ^. 

45x214 r +U[r ) 


In order to increase accuracy of estimate for E{r), we let p + g — 9/2 = 0, or 
q = 9/2 — p. Then we get 


/I / \ a n 7) (4p 11) 8 , ^ 10\ 

(1.15) -E{r) - Sg/ 2 -p,pil,r) = - ^ + O [r }. 

Further, taking p = 7/4 or 11/4 yields 

7 1 — 

(1-16) 5'ii/4,7/4 (I,?-') = 4., 7/4 := -^5 ( r ), 

and 

-F;(r) - 5 'ii/4 , 7 / 4 ( 1 ,/) = y X 221 ^^^ ^ ’ 

lim (-F;(r)-S'n /4 7 / 4 ( 1 ,/)) = « 0.00025614. 

r-S-l~ \7r / 77 11 

Letting p = 2 or 5/2 yields 

/4 1 _ p5/2' 2 

^ 5 / 2,2 (l,/)=(^j 


(1.17) 

and 


5 1-r'- 


:=A (/) 


^E{r)-S,/2,2{l,r') = + O ir^°) , 

(2 \ 2 16 

lim (-E(r) - 5*5/2 2 (1,/) ) =-— ft!-0.0033802. 

r^l- \7r ) 77 25 


These show that 5ii/4_7/4 (l,r') and 55 / 2,2 (Ij^O be excellent approxima¬ 
tions for the complete elliptic integrals of the second kind. The purpose of the 
paper is to prove this assertion. Our main results are contained in the following 
theorems. 
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Theorem 1. The function 




’ 11 / 4 , 7/4 ( 1 ) ’’0 

is strictly decreasing from (0,1) onto (11 (tt — 2) / (47r) , 1). Therefore, the double 
inequality 


(1.18) 


1 — /i + /iS'ii/4,7/4 (Ij f') < —E (r) < 1 — A + AS'ii/ 4 _ 7/4 (1, r') 


holds if and only if fJ, > 1 and A < 11 (tt — 2) / (47r) Ri 0.9993. 

In particular, we have 
(1.19) 

>5'll/4,7/4 < f-® (’’) < + '^^^^|^r^5'll/4,7/4 (1)’'0 < ||■'S'll/4,7/4 , 

where the coefficients 1 and 22/ (Jtt) r: 1.0004 are the best constants. 

Theorem 2. The function 

'' l-S6/2,2(l,r') 

is strictly increasing from (0,1) onto (1,25 (tt — 2) / (Ott)). Consequently, the double 
inequality 

(1-20) 1 — ^ + f.S^/2^2 (Ij'^’O < ^ I — r] + ?7<S'5/2,2 (1) ’’0 

holds if and only if f > 25 (tt — 2) / (Ott) r 1. 0094 and rj < 1. 

Particularly, it holds that 

(1-21) If >5'5/2,2 (1; r') < _2(8^^25) 25(w^ S'5/2,2 (1; ’’0 < (’’) < 5'5/2,2 (Ij 1^') j 

where the coefficients 25/ (Stt) r 0.99472 and 1 are the best constants. 


The paper is organized as follows. Some lemmas used to prove main results are 
presented in Section 2. The proof of Theorem 1 is complicated and longer, so it 
is independently arranged in Sections 3; while the proof of Theorem 2 is placed 
in Section 4. In Section 5, some interesting and applied corollaries involving the 
monotonicity of difference and ratio between ( 2 / 7 r) E (r) and Sg/ 2 -p^p (1, r') are de¬ 
duced. In the last section, it is shown that our approximations S' 11 / 4 ^ 7/4 (l,r0 
^ 2 , 5/2 (Ij’^O 8 ire indeed excellent by accuracy of comparing some known approxi¬ 
mations with ours. 


2. Lemmas 

In order to prove our main results, we need some lemmas. 

The first lemma is the ”L’Hospital Monotone Rule” (see [JS], [27]), which has 
been widely used very effectively in the study of some areas. A new and natural 
way to prove this class of rules can refer to [48] . which is easily understood and 
used. 

Lemma 1 ([46l Proposition 1. 1 ], m Theorem 2]). For —00 < a < b < 00 , 
let f,g : [a,b] —>■ ffi. 6 e continuous functions that are differentiable on {a,b), with 
f (a) = g (a) = 0 or f (b) = g (b) = 0. Assume that g'(x) 0 for each x in (a, b). 

If f /g' is increasing (decreasing) on (a,b) then so is f/g. 
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The second lemma is a monotonicity criterion for the ratio of power series, which 
will be used to prove Theorem 2. 

Lemma 2 (@1]). Let A (t) = J2T=o ^ (^) ~ power 

series converging on (—r,r) fr > OJ with bk > 0 for all k. If the seguence {ok/hk} is 
increasing (decreasing) for all k, then the function 1A (t) jB (t) is also increasing 
(decreasing) on (0,r). 


A more general monotonicity criterion for the ratio of power series has been 
established in recently. 

For a; > 0 the classical Euler’s gamma function F and psi (digamma) function if 
are dehned by 

r°° r' (t\ 

(2.1) r(a:)=/ t^-^e-^dt, = 

Jo r(x) 

respectively. In the proof of Theorem 1, we will use several inequalities for the 
gamma and psi functions, in which Lemma [3] is very crucial. 


Lemma 3 ([51], (52] (2.8)]). For all x > 0 and all a € (0,1), it holds that 

1 — a 


X 


X a 


Tjx + a) 
x^r{x) ’ 


or equivalently, 

( 2 . 2 ) 


F (a; + a) 

(x+ «)'-“ ^ r(x + i) 


< 


< 


f^l — a ' 


Lemma 4 ([S3]). For x € (0,1), it holds that 


+ 1 x^ + 2 

—— < F (x + 1) < X € (0,1), 

X 1 X -\-2 


or equivalently, 
(2.3) 


1 


< F (x) < 


x^ + 2 


,x e (0,1). 


(x + 1) X (x + 2) X 

Lemma 5 f [54l Lemma 1.7]). Let x be a positive real number. Then we have 


ifix + 1) > ln(x + -). 

Lemma 6 (jS^ Lemma 7], [IHl Lemma 1]). Let n G N and m G NU{0} with n > m 
and let Pn (t) be an n degrees polynomial defined by 


(2.4) Pn{t)= Y, 

i—0 

where a„, Om > 0, > 0 for 0 < i < n — 1 with i ^ m. Then there is a unique 
number tm+i G (0,c») to satisfy Pn (t) = 0 such that Pn (t) < 0 for t G (0,<m+i) 
and Pn (t) >0fortG {tm+i,oo). 
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3. Proof of Theorem [U 
We are in a position to prove Theorem [T] 


Proof of Theorem [IJ Let us consider the ratio 


1 - (2/77) E (r) 11 (1 - r'V^) (1 - (2/7r) E) /i (r) 

1 - 5 ii/4 .7/4 (1, r') 11 (1 - r'V4) -7(1- r'ii/4) ’ (r)' 


Clearly, /i (0+) = /s (0+) = 0. 
Differentiations by the formulas 


(3.1) 


dK _E-r' K dE _ E-K d{K - E) _rE 
dr rr'^ 'dr r ' dr r'^ 


yield 

f'Ar) ^ 1 - (2/t) E + 8 - r'^) {K - E) / {inr^) __ f, (r) 

/2 (r) 1-r' ' U (r) ’ 

/3(0+) = /4(0+)=0, 


777/3 (r) 

2 /I (r) 


j- I _ 7^2^ 7-'1/4 _ 8,,'2 _E — K 


p ( - r '2 

- 4-- E- 


y .12 


(K-E)- 7- 


(8 - 3r2) E- {8- 7r2) K 

j,4j,/3/4 


+ ((7r2 + 8) - (llr2 + 8) E) 


i!!! /3(o y 
2 


where 

/s (r) = 327s: - 327; - Idr^TsT - 3r^K - 2r^E, 

fe (r) = i (1287s: - 12871 - 224r27s: + 93r‘‘7s: + IGOr^Tl - 21r'‘7;) . 

If we can prove that fr (r) > 0 for r G (0,1), that is, the function f^/f^ is 
decreasing on (0,1), then so is f(/f 2 by LemmalU which in turn implies that / 1//2 
by using Lemma [T] again, that is, the function F is strictly decreasing on (0,1). 
Next we prove that fr (r) > 0 for r G (0,1) stepwise. 


32K - 32E - Ur^K - 3r*K - 2r^E 
^ 1 1287G - 1287; - 224r27s: + Odr-^Ts: + leOr^Tl - 21r'‘7; 
h (t-) - (t-) fi (r) 

r . ■ r . 1 
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Step 1: Expanding /r (r) in power series. By the expansions (|1.3p and (II.4p . we 
have 

/s (r) = i2K - 32E - lAKr^ - - 2r‘^E 

/ OO / 1\2 oo /'_1\ /I'l 

_ — I 32 \ ' ^ 2 / n ^2n ^2 \ ' ^ 2 J n \ 2 J n ^2n 

^ n!n! 


n—0 

oo 




n—O 


n\n\ 


J2n 


oo , 1.2 

_ O 2» 

^ om ' 


- 2 E 




^ ^ (n —l)!(n —1)! '^ / j (n—2)\{n—2)\' ^ / j (n —l)!(n —1)! 

n=l n=2 n—1 

^ ^ n(5n-6) (n - 1) (n - 2) 2 n ^ 6 ^ 2 n 

2 ^ n!n! ^ ' 2 ^ ^ ’ 

n=3 n=0 


where 

(3.2) 


(5» + 9)(i)"^i ^ 5n + 9 / r(n + 3/2) Y 
nl{n + 3)l nl{n + 3)l V r(l/2) J 


k k) = 


{123K - 128E - 224r2i4: + 93r‘‘i^ + IhOr^E - 21r^E) 

{-DM 

n—0 

2n 


/ 00 / 1 \ 

128 V 

n\n\ 


„2n 


n—0 


-128E 


/ n \ 'Z f n j,2n 


n\n\ 


-224 E 


( 1)2 

V2/n-l 


(n- l)!(n- 1)! 


+ 93E 


( 1)2 

\ 2 Jn -2 


2n 


n—2 

00 


(n-2)!(n-2)! 

{-MiD.- 


n—2 V2/n—2^2n 


(n-2)!(n-2)! 


°° (- 1 ) (i) 

+160 E , ^ ^ - 21 E 

1)1 in-1)1 

37 r “ n (n - 1) (n - 2) (n - 4) (n + 15) (2n - 5) 
~16 ^ 


j.lr?.! 


n-3^2n 


n—3 

Stt 


16 


'E< 

n—0 


2n 


where 

(3.3) 


_ (n-l)(n+18) (2n+l)(4)^ ^ (n - 1) (n + 18) (2n + 1) / r(n+l/2) y 
n! (n + 3)! n!(n + 3)! \ r(l/2) / 

Also, using binomial series we have 


(3.4) 


,,'-7/4 _ 


= (i-E“^'^ = E 


(I) 


n ^2n 


:=E^- 


n =0 


n—0 


Then making use of the Cauchy product gives 

2\-7/8 „ , . Stt 


OO oo oo 

n=0 n=0 n=0 


Stt 

16^ 


' E ( 8a„ + E 


, 2n 6 j 

fcCfc I r := —r 2_^ dn 


n—0 


fe=0 


16 


n—0 
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Step 2: A simple verification yields do = di = d 2 = ds = 0 and di, do, do, dr, 
ds, do, dio are equal to 

35 903 7343 193 225 36 001035 387471275 7897834 945 

32 768’ 262 144’ 1048 576’ 16 777 216’ 2147 483 648’ 17179 869 184’ 274 877 906 944’ 

respectively. 

Step 3: We prove that for n > 10, 


(3.5) Dn ■■= -dn+i - dn> g{n), 

where 


, X 4 n(2n+l) 3 128 , , 

( . ) gin) - - 2) (n + 3) ~ 7r (7/8) n + 1 77r 64n - 9^^ 

here 


(3.7) 


9i (^) = 

k=2 


n — k (fc — 1) (/c + 18) 
n — fc + 1 (fc + 1) (fc + 2) (fc + 3) 


n 

''y ^ Oln-kldk- 

k=2 


To this end, we note that for k >0, 


dfe+i 


bk+i 

Cfe+l 


1 (2fc + 3)^ i5k + 14) 
4(5fc + 9) (A: + l)(fc + 4) 
k + 7/8^ 

TTF'’*” 


5 


1 fc(2fc+l)(2fc + 3)(fc+19) ^ 
4(A:-l)(fc+l)(A: + 4)(fc+18) ^ ^ 


Also, it is seen that cq = —3, ci = 0, Cfe > 0 for k > 2, and Ofc, 5^ > 0 for /c > 0. 
Then the sequence Dn can be expressed as 


D 


n 


8 8 / \ f "" 

— y f 8c!.r),+ l H“ ^n+lCQ 4“ ^ ^ ^n+1 —fcCfc | I 80-7^ “h bfiCQ “h ^ ^ 

' \ fc -2 / \ k^2 



+ ^^OCn +1 + / ^ 
fc =2 


^^n+1 — bfi — k 


^k 


181 (2n + 3)^ (5n +14) 

1^7 4 (5n +9) (n + 1) (n + 4) 



3 


/8 n + 7/8 A 
V7 n+1 J 


b 


n 


8 1 n (2n + 1) (2n + 3) (n + 19) 

7 4 (n — 1) (n + 1) (n + 4) (n + 18) ^ 


/8n-fc + 7/8 
\7 n — k + 1 



bn — kC-k 


"Q-n 


8 n (5n^ — 6n — 29) 

7 (5n + 9) (n + 1) (n + 4)7 n 
2 n (2n + 1) (2n + 3) (n + 19) 

7 (n — 1) (n + 1) (n + 4) (n + 18) 


Cn + 


1 n — k 

7 n — fc + 1 

k=2 


bn—k^k' 


kCk 
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Notice that coefficient of a„ is positive due to bv? — 6n — 29 > 0 for n > 4, by 
the double inequality (12.211 in Lemma [3] it is derived that for k >0, 

1 (2fc + l)(5fc + 9) 

27r (fc + 1) (fc + 2) (/c + 3) ’ 

^ fci/8r(7/8)’ 

2 (fc-i)(fc + i8) 

^(fc + l)(A: + 2)(fc + 3) ^ - ’’ 

Applying these inequalities to the expression of yields 

^ ^ 8 n (5n2 - 6n - 29) 1 (2n + 1) (5n + 9) 3 u 1 

" ^ 7 (5n + 9) (n + 1) (n + 4) ^ (n + 1) (n + 2) (n + 3) ~ 7 u + 1 nVsr (7/8) 

2 n (2n + 1) (2n + 3) (n + 19) 2 (n — 1) (n + 18) 

7 (n — 1) (n + 1) (n + 4) (n + 18) tt (n + 1) (n + 2) (n + 3) 

("181 I ^ \ ^ ~ ^ _ ^ _2 (fc — 1) (fc + 18) \ 

7^ l^n-/c + l + 7/8)^/® TT (A; + 1) (/c + 2) (fc + 3)y ■ 

By the known inequality 

(1 — x)^ < 1 — px, p,x G ( 0 , 1 ) 


Ofe > 

1 

(fc + 7/8)^/® ^ 

Ck > 


it is deduced that for 2 < fc < n, 


1 


(n — fc + 7/8) 


1/8 ^ 


> 


(n - 2 + 7/8)^"'® ni/8(l-9/(8n))^/® 


= 64- 


,7/8 


(1 — 9/(64n)) 64n —9’ 

which is used to the last member of the right hand side in (13.81) and factoring gives 


Dn > 


n (2n -I-1) 


,7/8 


TT {n+l){n + 2){n + 3) 7T{7/8)n+l 


128 


77r 64n — 9 


—y — 

— Q \ n — 




n — k {k — 1) {k + 18) 


A: -I- 1 (A; -I- 1) (A: -I- 2) (fc -I- 3) / 


which proves the step. 

Step 4: The sequence gi {n) := ^n-kPk defined by (|3.7I) can be expressed 
as 


n(lln^+8n+2l) 


(3.9) 51 in) = i^ (n + 1) + 7) - ^^+ 4 ) ■ 

where ip (A) denotes the psi function, 7 is the Euler’s constant. 

In fact, decomposing rational function Un-kPk i^to partial fractions gives 

n — k (k — 1) {k + 18) 


C^n — kf^k — 


n — A: + 1 (A: -f 1) (A: + 2) (A: + 3) 


n + 1 1 

= -17-r X 


n -I- 2 A; + 1 

-30!i±5x ' 


n + 2 1 

48-T X 


Ti 3 k 2 
n{n + 19) 


-1-4 A;-I-3 (n-I-2) (n-I-3) (n-I-4) n — k + 1' 
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Hence, we get 


n I 1 / 1 


fc =2 


n + 2 \ k n + 1 2 

\fc=i / 


+48 


n + 2 


n + 3 \ k 

\k=l 


f n ^ 

Et 


1 


1 


-30 


n + 3 


/ n ^ 

Ei 


n + 1 n + 2 
1 1 


11 

Y 


n + 4 \ k 

\k—l 


+ 


25 


n + 1 n + 2 n + 3 12 


n{n+ 19) 


_ fy i-i 

(n + 2) (n + 3) (n + 4) 

+ 7n2 - 12 r + 24 ^1 
(n + 2) (n + 3) (n + 4) k {n + 1) {n + 2) (n + 3) (n + 4) ’ 


n (lln^ + 8n + 2l) 


which by the identity X]fe=i 'I = V' (''^ + 1) + 7 proves the step. 

Step 5: We show that the sequence g (n) defined by (13.61) satisfies the inequality 


(3.10) 


, , 128 n’’/® n® + 7n^ - 12n + 24 

Q (n) > ------- -02 (n). 

Ttt 64n — 9 (n + 2) (n + 3) (n + 4) 


for n > 4, where 


52 (n) 


In (n + 1/2) + 7 


n (lln^ + 8n + 2l) 

(n + 1) (n® + Tv? — 12n + 24) 


^ 8/— (n+4)(2n+l)(64n—9) (64n—9){n+2)(n+3)(n+4) 

22 (^+l)(^^+7n2 —12n+24) 128+(7/8) (^+r)C^^+'^^^“r2n+24) 


The inequality (I3.10p follows by the Lemma [SJ 
Step 6: We show that g 2 {x) > 0 for cc > 10. 

Differentiation yields 

! 7 (896a:'^+7346£c®+41033a;'^-3438£c^-149813£c3-227064£c^-40824a:+864) 

v^/ 256 a;7/8(a,4_[_8a,3_5a;2_|_i2a,_|_24)^ 

371 (55a:®+3806a;'^ + 17 10l£c^+216a:3-71 514x^-73 824a:-33 84o) 

128+ (7/8) (a:^+8a;3-5a:2 + 12a;+24)^ 

I ^ I lla:^ + 16a;^ + 182a:^ + 723;^-993a:^-384a;-504 
”^2a7 + 1 (a:<‘+8a:3-5£c2 + 12a:+24)^ ' 


Application of the second inequality of (12.31) in Lemma |4] gives 


r - < 


a;^ + 2 1 _ 177 6 

(x + 2) x_ 161 5 
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and since tt > 3, it is acquired that 
92 {x) > 


7 (896a;'^+7346a;®+41 033a:'^-3438£!;^-149 813£!;^-227 064£!;^-40 824a:+864) 

2^0 x'^/^{x'^+8x^—5x^ + 12x+24:)‘^ 

3x3 (55a;®+3806a;^ + 17101a;^+216a;^-71514a;^-738242;-33840) 

128 X 6/5 (x*+8x^ — 5x^ + 12x+24)''‘ 

^ _L llx^ + 16a;^ + 182a;*+72x^-993a:^-3843:-504 
(a:'‘+8a;3-5a:2 + 12a:+24)^ 


‘2x 1 


( Qiix) _ 


256(a;<‘+8a;3-5x2 + 12a;+24)^ ( 2 a;+l) 


TW ) > 


where 


53 {x) = 896x^ + 7346j;® + 41033a;'^ - 3438j;‘‘ 

-149 813x3 - 227 064x^ - 40 824x + 864, 

54 ( 2 ;) = 512x® + 15 474x^ + 153 661x® + 638 728x3 +482 131x‘^ 

-2496 996x3 _ 3737 3932-2 _ 2134 qqOx - 341712. 
Making a change of variable t = x — 3 > 7 yields 

53 (a;) = 896t^ + 26 162t6 + 342 605^® + 2450 487^"* + 10 008 901t3 

+22 815 555^2 + 26 081658t + 10 797192 > 0, 

54 (x) = 512^8 + 27 7Q2f + 607 639^® + 7103 356^® + 48 333 256^^ 

+ 194 587122^3 + 448 344153t2 + 530 387 220t + 235 084 068 > 0. 
Thus, to prove that 52 (x) > 0 for x > 10, it suffices to prove that 


35 (x) :=In 

Differentiation again leads to 

54 ( 2 ;) 


94 jx) 
2x -|- 1 


-In 7 


,93 (2;) 

rVS 


> 0 . 


95 ( 2 ^) = 


93 (x) , 7 1 


34 (x) 2x + 1 33 (x) 8 X 

1 x"^ + 8x3 _ ^^2 _|_ ^2x + 24 36 (x) 


X (2x + 1) 


93 (x) 34 (x) ’ 


where 


36 (x) = 6422 528x32 + 40 606 208x33 - 29 604 936x30 - 195 118 044xO 
-8157468 886x8- 54 727 744 833x3'- 28 074 816 632x0 
-33 746 602 635x0 - 132 036 870 576x4 _ yg 353 742 4743.3 


-16 962 249 960x2 _ 7592 953 7303 . _ gg ^ 424 . 

Lemma [5] implies that the polynomial qq (x) has a unique zero point xq € (0, 00 ) 
such that 36 (x) < 0 for x S (0,xo) and ge{x) > 0 for x G (xo,oo). This in 
combination with (7) = 56 640 373 211408 308 > 0 indicates that g^ (x) > 0 for 
X > 7. Therefore, 35 (x) > 0 for x > 7, and so 

1 
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which reveals that g 2 (x) > 0 for a: > 10. It follows that 

g 2 (x) > 92 (10) =7 + lnf - fit 0.037141 > 0. 

Step 7: Steps 6 and 5 show that <? (n) > 0 for n > 10, which in conjunction 
with Step 3 yield 

g 

Dn — ydn-i-i dji > 9 {ji) > 0, 

that is, dn+i > (7/8) dn for n > 10. Taking into account Step 2, we conclude 
that dn = 0 for n = 0,1,2,3 and > 0 for n > 4, and therefore, /y (r) = 
/s (u) — (r) > 0 for r S (0,1). Thus the function / 1 // 2 , that is, F, is 

strictly decreasing on ( 0 , 1 ). 

It follows from the decreasing property of the function F on (0,1) that 
—^ = lim F (r) < F (r) < lim F(r) = I, 

47r r—r— 


which implies inequalities (jl.lSp and the first and second ones in (jl.l9p . The third 
one in (|1.19p is eqnivalent to 


22-777 

477 


11(77-2) 

477 


5'ii/4,7/4 '5'i1/4,7/4 (1)^0 


ll(22-77r) 

287r 


(5'ii/4,7/4 (1) r') 


7_ 

11 


) < 0 , 


where the inequality holds due to the increasing property of Stolarsky means in 
their variables (PI) and r' G ( 0 , 1 ). 

Since 

,. (2/77) £;(r) ^ (2/77) i?(r) 22 

r^0+ Sii/4,7/4 (I,r') r^l- 5 ii/ 4_7/4 (I,r') 777 


the coefficients I and 22/ (Ttt) are the best. 
Thus we complete the proof. 


□ 


Remark 2. From the Step 7 in previous proof, we have dn+i > {7/8)dn for n > 
10. This is also valid for A < n < 9 by an easy verification. It is derived that 
dn > (7/8)” di for n> 5, and so we get 


n—0 n—4 



35 r® 
4096 8 - 7r2 ■ 


Thus we obtain that 


h {r) = h (»’) - r' (r) 


16 ^ ” 216 8-7r2- 

n=0 


4. Proof of Theorem [2] 

Now we prove Theorems 
Proof of Theorem\M Utilizing (jl.4ll gives 


2 . ^ 
1- E{r) 

TT 


E- 

n—l 


nln\ 


ri ^2n 
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Applying binomial series to (jl.l7|) we have 

'4 1 _ /5/2 \ _ 16 (1 - - 2 (1 - + 1 


5'5/2,2 — 


5 l-r'2 
16 1 ' “ 


25 


oo (_ 5 \ OO (_ 5 \ N 

E V 2/n ^2n _ 2 \ ' V i/n ^2n _|_ ^ 


1 --E 

25 ^ 


n—0 n =0 


n—1 


(n + 2 )! 


"+2 ^2 ti 

' 1 


which implies that 


1 - Ss,,,, (i.r') = ^ f; 


25^ (n + 2 )! 

n—1 


-E 




Thus, the function G (r) can be expressed as 

l-( 2 / 7 r)ii;(r) 


G(r) = 


E OO 

n=l 


1-55/2,2(1,?'') Er=iWnr 


2n ’ 


where, by the formula (a)^ = a (a + Vn and u„ can be written as 

6 1 /1\ 2 1 /3\ 

“ 5 (n + 2)! UE-E5(n + 2)! UJn' 


By Lemma [2J due to n„ > 0 for n > 1, to prove the function G is increasing on 
(0,1), it suffices to prove the sequence {u„/m„} is increasing for n > 1. In view of 
Un > 0 for n > 1 , which suffices to check that (vn+i/vn) Un — Un +1 > 0 . 

A simple verification yields 

Vn+i _ (n - 1/2) (n + 1/2) 

Vn (n + 1)^ 

and 

Vn+i 3 (3n + 1) /1\ 1 n^ —lln —6 /3\ 

Vn ^ 5 (n + 1)^ (n + 3)! \2/„ 10 (n + 1)^ (n + 3)! V^/„ 

Direct computations give 

t’n+i ,, _ = n n 3 21 47 1881 157531 

un+1 81920’ 512 000’ 1310 720’ 64 225 280’ 6710 886 400’ 


42 559 507 577 997 177 20 743 573 

2264 924 160 ’ 33 554 432 000 ’ 81 201 725 440 ’ 2061 584 302 080 

for n = 1 , 2 , ..., 11 , respectively. And, for n > 12 , it is evident that {vn+i/vn) n„ — 
Un+i > 0 due to n^ — lln — 6 > 0 . 

Consequently, we obtain 

1 = lim G (r) < G (r) < lim G(r) = 

r—>0+ r—^1“ QtT 


which implies inequalities (11.201) and the first and second ones in (ll.2ip . The third 
one in ( 11 . 211 ) is equivalent to 

2 ( 87 r- 25 ) I 25(7r-2) o /, /■> 25 o /l _ 25 ( 87 r- 25 ) fn /i /\ 16\ ^ n 

- ^ - 9 ;;—*2,5/2 (^,r j-^*2,5/2 U,r j - - ^ - (*2.5/2 > U, 
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where the inequality holds due to the increasing property of Stolarsky means in 
their variables (PI) and r' G (0,1). 

Since 


lim 


{2/tt) E(r) 
0+ Sb/2,2 (1) f') 


= 1 and 


lim 

r—^l~ 


{2/TT)E{r) 

^5/2,2 (l,u') 


25 

Stt ’ 


the coefficients 1 and 25/ (Stt) are the best. 
This completes the proof. 


□ 


Remark 3. Denote by Wn = Vn — Un- It is easy to verify that the relation 

Vn+l n„+i 

Wn — Un Un-^1 

Vn Vn 

holds for n > 0. From the proof of Theorem\^ we clearly see that 
(4.1) > 0 

Vn 

for n > 3, which, due to wi = W 2 = W 3 = 0, W 4 = 3 x 2“^"‘/5 > 0, means that 
Wn > 0 for n > 4. This also yields 

^ 00 00 CX3 

'S'5/2.2 ( 1 , r') - -E (r) = ^ ^ X! 

n —1 n —1 71—4 

so we have 


Proposition 1. The function 


Gi (r) = 


85/2,2 {iy)-iEir) 


is convex and strictly increasing from (0,1) onto ((3/5)2 16/25 — 2 / 7 r). Con¬ 

sequently, we have 

3 S 5/2 2 (I,r')-( 2 /TT)E(r) 16 2 

- X 2-^^ < ^ ^ ^ ’ — y < - a; 0.003802. 

5 r® 25 TT 

Remark 4. Further, the relation ra also indicates that for n > 4 

W4 3 X 2-14/5 3 

Wn > - Vn = - rrVn = - Vn , 

V4 175 X 2-14 ^ 875 

and therefore. 


9 °° O OP 

^ 5 / 2.2 ( 1 , r') --G(r) = 


Vnr^"^ = - 

, 875 875 , , 

—4 n—4 \n—l 


E - E 


VnV 


2n 


3/1 2 ^ , /I 2 84 5 e 

= r +«’' + 2 S’' 


which can be stated as a proposition as follows: 
Proposition 2. For r G (0,1), it holds that 

2 „ . . 875 „ 3 1 


’ cn / X 875 no 3 / 


, * 2 84 5 g 

1 - r - r-r 

4 64 256 
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5. Corollaries 

As direct consequences of Theorems [T] and [21 we have 


Corollary 1. Both the functions 

o . , iV^)E{r) 

r hA-A (r) - * 11 / 4 , 7/4 (l,r ) and r ^ --—- 

TT *11/4,7/4(1,?') 

are strictly increasing from (0,1) onto (0, 2 / 7 r — 7/11) and (1, 22/(Ttt)), respee- 
tively. And therefore, we have 

(5.1) 0 < -i; (r)-* 11/4 7/4 (1,/)<--;^R:!0.00025614, 

TT ' TT 11 


(5.2) 


1 < 


(2l-K)E{r) 22 

^ ^ ’ < — « 1.0004. 


5'ii/4,7/4 (I,?"') 771 


Proof. By the increasing property of Stolarsky means in their variables and r' = 
\/\ — r'^, it is seen that *11/4,7/4 (IjfO is decreasing with respect to r on (0,1), and 
so both the functions 

r ^ (1 - * 11 / 4 , 7/4 (1, r')) and r ^ ^ , 

* 11 / 4 , 7/4 U, ^ ) 

are positive and increasing on (0,1). 

From Theorem[l] we see that r i-A (1 — F (r)) is positive and strictly increasing 
on (0,1). It follows from the identity 

If (r) - * 11 / 4 , 7/4 (1, r') = {1-F (r)) (l - * 11 / 4 , 7/4 (1, r')) 

that the function r ( 2 / 77 ) F(r) — * 11 / 4 , 7/4 ( 1 , r') is also positive and strictly 
increasing on (0,1). Therefore, the double inequality (15.11) is valid. 

Making use of the assertion proved previously, and noting that 


( 2 / 77 ) F(r) 
* 11 / 4 , 7/4 (1, t') 


^ ~ * 11 / 4 , 7/4 (1, 


*11/4,7/4 (1, T'O 


gives the desired assertion. Then the estimate inequalities (j5.2l) follow. 


□ 


Using the same technique we can prove 


Corollary 2. Both the functions 


, , ( 2 / 77 )F(r) 

r -F(r) - * 5 / 2,2 (l,r ) and r ^ --—- 

77 *5/2,2 (1,7' ) 


are strictly decreasing from (0,1) onto (2/77 — 16/25,0) and ( 25 /( 877 ),!). 
therefore, we have 


(5.3) 

(5.4) 


- 0.0033802 
0.99472 


2 16 2 

---<-F(r)- *2,5/2 (l,r')<0, 
TT 25 TT 

25 ( 2 / 77 ) F(r) 

877 *5/2,2 (1,7'') 


And 


Now we give a monotonicity property for the ratio {2/tv) E{r) /* 9 / 2 -p,p (1, ?'') 
with respect to r for p G (— 00 , 9/4]. To this end, we need a known statement 
proved in [4TJ Theorem 5] by Losonczi (see also [56l Theorem 3.4]). 
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Lemma 7. For fixed c>0,0<x<y<z, the function 

‘^ 2 c— y) 


p I—>■ 


•— di2c—p,p Vi 


S2o-p,p(x,z) ■ 

is strictly decreasing on (—oo,c] and strictly increasing on [c, oo). 

Corollary 3. Let p G (— 00 ,9/4]. The the function 

{2/TT)E{r) 

^ ^ 

*^9/2—p,p ) 

is strictly increasing from (0,1) onto (1,2/ (TrOp)) if and only if p G (— 00 , 7/4] and 
strictly decreasing from (0,1) onto {2/ (nOp) ,1) if p G ]2,9/4]. Consequently, we 


have 


(5.5) 

*^9/2—p,p (I 5 ^ ) 

(5.6) 

■^Sg/2-p,p{iy) 


TT 

2 

—1 

TT 


2 

ttO 


p 


where the coefficients 1 and 2/ {nOp) are the best possible, here Op is defined by 

prm - 

Proof, (i) The necessity can be derived from 

lim > 0 . 

r->-o+ 8r' 

From limr_^o+ Rp (^) = liinr^o+ Sg/ 2 -p,p (1,t') = 1 and L’Hospitial rule it is ob¬ 
tained that 


1 = lim 

= lim 


In Rp (r) 


1 


r^ 0 + \Rp (r) - 1 Sg/ 2 -p,p (1, r') 
In Rp (r) 


= lim 


d (lni?p (r)) /dr 


^ 0 + {2/'n:)E{r) - Sg/ 2 -p,p (l,r') 7 -^ 0 + d ((2/7r) £'(r) - Sg/ 2 -p,p (l,c0) /dr 

By the expansion (11.151) and L’Hospitial rule we have 

(4p - 7) (4p - 11) {2/Tr)E{r) - Sg/ 2 -p,p{l,r') 

- 71 - = iim -- 

5 X 214 ^^ 0 + p8 


= lim 

r—>-0+ 


d {{2 /tt) E (r) - Sg/ 2 -p,p (1, r')) /dr 


It follows that 

d {\nRp (r)) /dr 


lim 

r—>-0+ 


Srl" 


= lim 


d (In i?p (r)) /dr 


d (( 2 / 77 ) E (r) - Sg/ 2 -p^p (1, r')) /dr'' 


0 + yd{{2/TT)E{r) - Sg/ 2 -p,p{l,r')) /dr 
(4p-7) (4p-ll) 


8r7 


5 X 214 

which together with p G (— 00 , 9/4] gives the necessary condition p G (— 00 , 7/4]. 

(ii) To prove that the condition p G (— 00 , 7/4] is necessary, we have to prove 
that the function 

'S' 9 / 2 —po.PO l" ) . 7-)* /.J\ 
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is strictly increasing (decreasing) on ( 0 , 1 ) if p < {>)po- Dae to r' = \/l — it 
suffices to prove that r' i- 7 > Rpg,p{f') is strictly decreasing (increasing) on ( 0 , 1 ) if 
p < {>)pq. Assume that r[, G (0,1) with r[ < Then 1 < l/r^ < l/r). By 
Lemma [71 we have 

S 2 c-po,po ( 1 , 1 A^) . . ^ S 2 c-p,p{IMt'^) 

S2c-p,,pA^Mr[) ^ ’ S2c-p,p{lMT{y 
which is equivalent to 


S2c—po,po (1,1/ ^ 2 ) / \ ^2c-po,po (l,l/^l) 

52c-p.p(l,l/r') ^ ^ ^ 2 c-p,p ( 1 , 1/r'i) • 


From the homogeneity and symmetry of Stolarsky means with respect to their 
variables, this shows that r' !->• Rp^^p (r') is strictly decreasing (increasing) on ( 0 , 1 ) 
if p < (>)po- 

Now, if po = 7/4 and p G (— 00 , 7/4), then r !->• Rp^^p {r') is positive and strictly 
increasing on (0,1). While Corollary |T] tells us that r !->• Rpo (r) is also. It is 
deduced by the relation 


Rp{r) 


{2/TT)E{r) ^ Sg/ 2 -pa,po (l.'^O 

5'9/2-po,po (l,r') , 59 / 2 -p.p(l,r') 


:= Rp, (r) X (r') 


that r ^ Rp (r) is also strictly increasing on ( 0 , 1 ). 

(iii) We continue to show that the function r ^ Rp (r) is strictly decreasing if 
p G [2,9/4]. Let po = 2 and p G (2,9/4]. Then r !->• Rpg^pir') is positive and 
strictly decreasing on (0,1). Similarly, this together with Corollary [7] reveals that 
r ^ Rp (r) is also strictly decreasing on (0,1). 

(iv) Lastly, the inequalities (15.51) and (15.61) follow from the monotonicity of the 
function r ^ Rp (r) on ( 0 , 1 ). 

The proof is finished. □ 


Remarks. Let p = 9 j A; 2^9 j A in Corollary\^ Then hy the monotonicity 

of S2c-p,p and y/Op) S 2 c-p,p in P on (0,c) given in (P2) and (P3) we have 

7769/4 (l,r') < A 3/2 (l,r') < 5 'ii/ 4 y /4 (l,r') < ;^A(r) 

22 , 2^/3 , 2 X 3"^/® 

< :;-5'ii/4,7/4 (1, ) < - A3/2 (1, r ) < - Hcg/i (1, r ), 

tTT TT TT 

2p4/9 25 2 

—~—7g/4 (1, r') < — S' 5 / 2,2 (1, t') < —R (^) < ^5/2,2 (1, 1"') < Ig /4 (1, r') , 

where Hcp (a, h) = Pie (a^, and Ip (a, b) = I (a^, are the p-order Hero- 

nian mean and identrie (exponential) mean of positive numbers a and b, respectively. 

Using expansion (11.151) and Corollaries IH together with the property of Sto¬ 
larsky means (P2), we obtain immediately 

Corollary 4. For p G (— 00 ,9/4], the inequality 
(5.7) ‘^E (r) > Sg/ 2 -p,p (1, r') 

holds for all r G (0,1) if and only p G (— 00 , 7/4]. The inequality reverses for 

pG [2,9/4]. 
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Remark 6. Taking p = —9/2,—9/4,0, 9/8,3/2, 7/4; 2,9/4 in Corollary^ we get 
immediately 

(1, rO Heg/2 ( 1 , r') < L9/2 (1, /) < Heg/4 (1, /) < 

(5.8) ^3/2(1,'^') < Sii/4^r/4{l,r’) <^E{r) < S5/2.2i^,r') < Ig/4{l,r') 

holds for r G (0,1), where Lp (a, h) = L (qP, is the p-order logarithmic mean 

of positive numbers a and b. 

Remark 7. For a,b > 0 with a ^b, the Toader mean T(a, b) is defined in [ST] by 

2 , - 

T{a,b) = — / V cos^ t + 6^ sin^ tdi. 

TT Jo 

An easy transformation yields 

\ {b/af\ ifa>b, 

T{a,b) = \ ) - { 

y He [yl- (a/b) j ifa<b. 

Thus, all our results can be rewritten in the form of Toader mean, for example, 
inequalities \ 1.1 91 \ 1 . 21 \ and \ 5 . 8 \ are equivalent to 

5'ii/4 ,7/4 {a, b) <T (a, b) < 22^ max (a, &)+^^^^/^S'ii/4,7/4 (a, b) < ||S'ii/4,7/4 (a, 
H 5'2,5/2 (a, b) < max (a, b)+ '^^^g~'^'> 82^5/2 (a, b) <T (a, b) < 82,5/2 (a, b ), 


^9/2 (®>(o, 5 ) < Y'G'(a^' 6 )iie^^^ 7 (a^ < ^9/2 (a,&) < ^69/4(0,6) < 

^3/2(0, &) < -S'11/4,7/4 (a,^) < T(a, &) < 5 ' 5 / 2,2 (a, &) < T9/4 (a, &) ■ 


6. Comparisons with some known approximations 

Let .4(r) be the given approximation for ( 2 / 77 ) ^^(r) and let A (r) = A{r) — 
( 2 / 77 ) E (r) denote the error. In general, the most important criterion to measure 
of accuracy of the given approximation A (r) for ( 2 / 77 ) E (r) should be the maximum 
absolute error maXrg(o,i) |A (r) |. Due to r G (0,1), however, similar to Barnard et 
al.’s opinion in [9], if A (r) = + O (^A'^0+2'^ with Eng 7 ^ 0 by expanding in 

Maclaurin series, then the leading item i5o := £noT^"'° can be viewed as a measure of 
accuracy of the given approximation A (r) for ( 2 / 77 ) E (r). For this reason, we call 
A{r) an ng-order approximation for {2/Tr) E{r). And, .4(r) is called an ng-order 
lower (upper) approximation for ( 2 / 77 ) if (r) if A (r) < (>)0 for all r G (0,1). In 
most cases, the greater the order rin is, the higher the accuracy of approximation 
A (r) for ( 2 / 77 ) E (r) is. 

Of course, a desirable approximation A (r) also has the simplicity of the expres¬ 
sion. Unfortunately, there exists frequently, certain negative correlation between 
the accuracy and simplicity. 

Now we choose those approximations which have higher accuracy mentioned in 
Introduction to compare with our ones. 
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6.1. For some lower approximations. In Tables 1, the values in the third col¬ 
umn are derived by expanding in power series, while max,re(o,i) l^i (?’) I for z = 1, 
2, 5 in the fourth column are from [8l Theorem 1. 1], [211 Theorem 2.5], (15. ip . 
respectively. 


Table 1: The lower approximations 


{r') 

Expressions 

C/Q •— tno ' 

maXrg(o,i) |A* (r) ] 

Ml (r') 

^ 3/2 ( 1 ) A) 

1 ^8 

= 2 _ 2-2/3 ft! 0.0066592 

TT 

^2 (r') 

23A(iy )-sff(iy )- 2 s(iy) 

3 „12 

2 2 S- 2 V 2 ^ r. nofiofisi 

16 


-43 ir') 

J (9r'^ + 14r'+9)^ 

128 pTI? 

1 „12 

> 0.0038073 

A 4 (r') 

l+r'+r'^ 1 1+r' 

2(1+9-') 8 

263 „8 

jT^r 

> 1 - 1 Ri 0.011620 

— TT 0 

As (A) 

^11/4,7/4 (1) A) 

I zn 

7x221 ' 

= ^ « 0.00025614 

77 i i 


Moreover, we can prove 


Lemma 8. Let Ai {x) (i = 1, 2, 3,4, 5j fee defined on (0, oo) by U.5]} . U.10\) . 
il.lAl and respectively (see also Table 1). Then the inequalities 

(6.1) As (x) > As (x) > A 2 (x) > Ai (x) > A 4 (x) 

hold for X > 0 with x 1. 

Proof, (i) The first inequality in (16.11) is equivalent to 


11 (a;’’-1) 128(x4-hl)^ 


Factoring yields 


As {x'^) - As (a:^) = 


1 


(a: - 1) 


rhi (x) > 0, 


1408 — 1) {x'^ + 1)^ 

where 

hi (x) = 5x1® 35^15 420 x 1® -f 966x1^ -f 1722xii -f 2268xi° -f 2415x® 

-f2362x® -f 2415x1' -f 2268x® -f 1722x® -f 966x^ -f 420x® -f 140x^ -f 35x -f 5. 
(ii) A direct computation leads to the second inequality in (16.11) is equivalent to 

V 2 


As (x) - As (x) = 


23^-5l^-V^ 


1 (9x^-I-14x-I- 9 )^ 
m (x + i)® 16 

(1 -I- x^) 1 llx^-I-36x®-I-34x^-I-36x-I-11 


16 

where the inequality holds due to 

2 


128 


(x-f 1) 


> 0 , 


/ a/2 (1 -f x^) \ / 1 36x-I-34x^-I-36x®-I-llx"i-I-11 

16 ) ~ (x + 1 )® 

1 (x - 1 )® (7x2 332 . 7 ) 


> 0 . 


(x + l)° 

(hi) It has been shown in |2U Lemma 3.2] that the third inequality in (16.11) holds 
for X > 0 with x / 1 . 
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(iv) The last inequality in (16.11) is equivalent to 



{x- 1)® (3a;6 + 18x5 - + 28x3 - Sx^ + 18x + 3) 

512 (x 2 _|_ x)3 

This completes the proof. □ 

Remark 8. In Table 1, both Ai {r') and A 4 {r') are A-order lower approximations. 
And for the accuracy, since < £ 4 ^^ and maXrg(o,i) |Ai (r) | < maXrg(Q 4 ) IA 4 (r) |, 
so the former is better than the latter. And, the last inequality in 16.11) also proves 
this assertion. 

While A 2 (r'), A^ (r') and A 5 {r')(= <S'ii /4 7/4 (l,r')J are 6 -order lower approxi¬ 
mations. In view of < £g^^ and 

max lAg (r) I <C min ( max IA 2 (r) I, max lAg (r) I I , 

>•6(0,1) \rG(0,l) >•6(0,1) / 

we claim that the accuracy of A 5 {r') (= i 5 'ii/ 4 , 7/4 {l,r')) is far superior to A 2 (r') 
and A 3 {r'). And, the first and second inequalities in a confirm similarly prove 
this claim. 

To sum up, our lower approximation A 5 {r')(= S'11/4,7/4 {l,r')) for { 2 / tt ) E (r) 
is the best of all five ones listed in Table 1. 

6.2. For some upper approximations. In Tables 2, the values in the third col¬ 
umn are derived by expanding in power series, while max,re(o,i) |Ai (r) | for i = 7, 
8 in the fourth column are from [SHI Theorem 3.2], (|5.3I1 . respectively. 


Table 2: The upper approximations 


A, (r') 

Expressions 

p(^)^2no 

max^g(o.i) |A(r) | 

A (r') 

>Cl/4 (l,u') 

J- 

> 1 - ^ Ri 0.36338 

A(r') 

18A(l,r')-5G(l,r')+3S(l,r') 

7 „12 
-^r 

18+3\/2 2 Q 058461 

16 


A (r') 

>5'5/2,2 (1, r') 

5x210 ' 

= ^ A 2 0.0033802 

25 TT 


Further we have 


Lemma 9. Let Ai (x) (i = 6,7,8) be defined on (0, 00 ) by S1.8\) . 11..91) . ([1.17| ), 
respectively (see also Table 2). Then (i) the inequality 

(6.2) max {A 7 (x), As (x)) < Ae (x) 

hold for X > 0 with x ^ 1 ; (ii) there is a xq € ( 0 , 1 ) such that .Ag (x) > A 7 (x) for 
X € ( 0 , xo) and As (x) < A 7 (x) for x € (xq, 1). 

Proof, (i) It has been proven in uni Lemma 5.1] that .A 7 (x) < Aq (x). To prove 
.Ag (x) < Ae (x), it suffices to prove that .Ag (x^) — As {x'^) < 0. Straightforward 
calculation gives 

A (A) - A (A) = 16(x3°-1)" _ x5 + l ^ (x - 1)® (x + 1) (x5 + 1) 

^ ’ 25(x8-l)" * + l 25(x8-l)' 

X (9x® -f 20x^ -f 44x® -f eOx® -f 74x^ -f GOx® -f 44x2 -f 20x -f 9) , 
which is obviously negative for x > 0 with x ^ 1 . 
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(ii) We have 

As (x^) - At (x^) = 


16(x^-iy foi + x 


2 5 

— —X - 


3 1-hx^ 


(6.3) 
where 
/i 2 (x) = 


25(x‘^-lf [8 2 16 ' 16 

/i 2 (x) — hs (x) 1 /i| (x) — hi {x) 


400 400 /i 2 (x) + hs (x) ’ 

31x® + 187x^ + 118x® + 49x5 + 430x4 + 49x3 + llSx^ + 187x + 31 


(x^ + 1 )^ (x + 1 )^ 


hs (x) = 75 


1 +X4 


Making a change of variable u = x + 1/x yields 
31^4 + 187^3 - 6^2 - 512u + 256 


/l 2 (x) = X 
and therefore, 
(6.4) 
where 


{u + 2 ) 

h'l (x) - hj (x) = - 


and hs (x) = 75x 


m2 - 2 


2/.„n 1,2 / a;2 { u - 2 f 


2 {u + 2 )' 


■hi (u), 


hi (u) = 3703m5 + 14 124m5 - 16 260^4 - 98 560m3 - 23 040m^ + 98 304u - 32 768. 
Since m = x + 1/x> 2, replacing m by x + 2 leads us to 

hi (x + 2) = 3703x®+58 560x5+347160x4+928 800x3+1014 000x^+144 OOOx-288 000, 


where x > 0 . 

It follows from Lemma [ 6 ] that the polynomial hi (x + 2) has a unique zero point 
xi G ( 0 ,oo) such that ^4 (x + 2 ) < 0 for x G ( 0 ,xi) and /14 (x + 2 ) > 0 for 
X G (xi,oo). Numeric computation gives xi G (0.399475162,0.399475163). This 
together with (16.41) and (16.31) reveals that 



for X G (0, xi), 
for X G (xi, 1), 


where Xi + 1/xi = 2 + xi, that is, Xi = ^xi + 2 — y^xi (xi + 4)j /2 

which proves the second assertion, where xq = xf ~ 0.28825. 

This completes the proof. 


0.53689, 

□ 


Remark 9. From Table 2, as upper approximations, A 6 (r'){= Ci/i{l,r')) and 

A 8 {r')(= Ss/ 2,2 have 8 -order accuracy, but the facts |£ 4 *^| < and 

maXrg(Q 4 ) jAg (r) | < maXj.g(Q 4 ) jAg (r) | together with the inequality i 6 . 2 \) show 
that the accuracy of As {r') is higher than As (r'). 

The upper approximation A 7 (r') has 6 -order accuracy, but its maximum absolute 
error is greater than our approximation Mg (r') ’s. By the second assertion of Lemma 
0 it is seen that there is a unique xq G (0,1) such that (2 /t:) E{r) < At (r') < 
As{r') for r G (0,ro) and {2/-K)E{r) < As{r') < At {r') for r G (xq, 1), where 
ro = v'l-x^ Ri 0.95756. 
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Remark 10. Not only that hut our approximation S' 11 / 4 ^ 7/4 (l,r^) and S 5 / 2.2 

have very small absolute relative errors. Exactly, from Hl.lfA) and \1.21\) it is derived 

that 


(2/7r) E (r) - Sn/ 4 , 7/4 (1, r') 


(2/tt) E Ir) 

(2/77) i?(r)-^ 5 / 2,2 (1,0 


< 1 - 


Ttt 


0.00040234, 


< 


Stt 


- 1 Ri 0.0053096. 


( 2 / 7 r) E (r) 

Furthermore, by Corollary the absolute relative error function 

(2/tt) E (r) - Sg/ 2 -p,p (1, r') 


£p (r) = 


{2/TT)E{r) 

is strictly increasing in r from (0,1) onto (0, |1 — 7r0p/2|) forp G (0, 7/4]U[5/2, 9/4]. 
Also, by the Remark{l\ the maximum absolute relative error inax^g(Q 4 ) fp (r) = 
|1 — 7r0p/2| is strictly decreasing in p on (0, 7/4] and strictly increasing on [2, 9/4]. 

Lastly, we close this paper by proposing a conjecture as follows. 

Conjecture 1. Let po ~ 1.763135 denote the unique solution of the equation Op = 
2/tt, that is, 

P J tt 

on (0,9/4]. Then there is a rg G (0,1) such that the function 

H{r)= 

1 — 5'9/2-po.po ( 1 , 

is strictly increasing on (0,ro) and strictly decreasing on (ro,l). Consequently, it 

2 

-E{r) < Sg/2-po,po (1,^') 


holds that 


for r G ( 0 , 1 ) with the best constant pq. 
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